Tunable multi-photon Rabi oscillations in an electronic spin system 
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We report on multi-photon Rabi oscillations and controlled tuning of a multi-level system at 
room temperature (S = 5/2 for Mn 2+ :MgO) in and out of a quasi-harmonic level configuration. 
The anisotropy is much smaller than the Zeeman splittings, such as the six level scheme shows 
only a small deviation from an equidistant diagram. This allows us to tune the spin dynamics by 
either compensating the cubic anisotropy with a precise static field orientation, or by microwave 
field intensity. Using the rotating frame approximation, the experiments are very well explained by 
both an analytical model and a generalized numerical model. The calculated multi-photon Rabi 
frequencies are in excellent agreement with the experimental data. 
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PACS numbers: 03.67.-a 71.70.Ch 75.10.Dg 76.30Da 

I. INTRODUCTION 

Harmonic systems are basic manifestation of quantum 
mechanics and appear in various forms, down to the 
nanoscale, as electro-magnetic or mechanical oscillators. 
In the case of a finite number of excited states, multi-level 
systems are proposed to perform quantum algorithms in 
size-limited 1,2 or scalable 3 schemes by using microwave 
(MW) pulses to generate entangled states. The coherent 
manipulation of spins in a multilevel system is fundamen- 
tal to the implementation of the Grover algorithm 4 . At 
the same time, the quasi-harmonic nature of the system 
can lead to interesting effects, where the two- or multi- 
level nature of a system can be interchanged and tuned 5 . 

Spin systems benefit from relatively large coherence 
and relaxation times, which makes them suitable as quan- 
tum bit implementations or as a new type of quan- 
tum random access memory 6-9 . In some studies, spin 
qubits operation is demonstrated at temperatures up 
to ambient 10 . Recent studies demonstrated that in di- 
luted systems, spin-spin dipolar interactions are suffi- 
ciently low to allow coherent, quantum manipulations. 
Such situations arc well exemplified by systems like 
the nitrogen-vacancy color centers in diamonds 11 ' 12 , N 



atoms in Cgo 13 , Er 3+14,15 and Cr 5+16 ions and molecular 
magnets 17,18 . 

In a previous work 5 , we demonstrated multi-photon 
spin coherent manipulation in a multi-level system 
(Mn 2+ ) diluted in MgO, a highly symmetric nonmag- 
netic matrix. In this work, we report on the possibil- 
ity of using a multi-level system to tune the nature of 
quantum Rabi oscillations, using a combination of two 
parameters: microwave power and/or magnetic field ori- 
entation. To this end, one needs a magnetic system with 
a well defined anisotropy, but sufficiently small in size 
to be overcome by microwave amplitudes achievable in 
typical experimental conditions. 



The S=5/2 Mn 2+ spin is our system of choice. As 
detailed in the next section, the crystalline anisotropy for 
this high cubic symmetry is orders of magnitude smaller 
than the magnetic (Zeeman) energy and therefore the 
multi-level system is quasi-harmonic. This is essential for 
successful multi-photon spin manipulation and for state 
tunability by magnetic field orientation. At the same 
time, the anisotropy remains smaller than or comparable 
to the microwave drive, an aspect that is essential for 
tunability by drive. 

In this S=5/2 system, we demonstrate multi-photon 
Rabi oscillations and controlled tuning of the system in 
and out of a quasi-harmonic level configuration. Using 
the rotating frame approximation, the experiments are 
very well explained by both an analytical model and 
a generalized numerical model. The calculated multi- 
photon Rabi frequencies and amplitudes are in excellent 
agreement with the experimental data. 

The article is structured as follows: in Sec. II we 
describe the spin Hamiltonian and its parameters; in 
Sec. Ill the experimental procedure and setup is detailed; 
in Sec. IV an analytical and numerical model is given, de- 
scribing the tunability of the multi-photon dynamics and 
in Sec. V we describe the experimental findings. 



II. QUASI-HARMONIC SPIN HAMILTONIAN 

The spins S—5/2 of the Mn 2+ ions are diluted in a 
MgO non- magnetic matrix of cubic symmetry F m ^ m (lat- 
tice constant 4.216 A). The Mn 2+ ions are located in 
substitutional positions of Mg 2+ ions. The high degree 
of symmetry ensures that the spins are seeing an almost 
isotropic crystalline environment and thus the fourth or- 
der magnetic anisotropy can be made much smaller than 
the Zeeman splittings. Interactions between Mn 2+ ions 
and symmetry deformations are neglected. The spin 



Hamiltonian at resonance is given by 19,20 : 

H = a/6 [Si + S* + St- S(S + 1)(3S 2 - l)/5] (1) 
+~fH Q ■ S-AS-I + ~fh mw ■ S cos(2ir ft) 

where 7 = gus/h is the gyromagnctic ratio (g = 2.0014 
the g-factor, /Xb— Bohr's magneton and h— Planck's 
constant), S XlVlZ are the spin projection operators, S 
is the total spin, a = 55.7 MHz is the anisotropy con- 
stant, A = 244 MHz is the hyperfine constant of 55 Mn 
(/ = 5/2), h mw and / represent the MW amplitude 
and frequency respectively, and Hq is the static field 
(Ho -L h mw ). In our experiments, the static field ensures 
a Zeeman splitting of 7-ffo ~ / ~9 Ghz, much stronger 
than all other interactions of Eq.(l). This implies that (i) 
.ffo's direction can be approximated as the quantization 
axis and (ii) coherent MW driving is confined between 
levels of same nuclear spin projection m/ (see also refs. 
5, 21, and 22). 

The six-level system is consequently quasi-harmonic, 
as shown by the level diagram in Fig. 1. The dashed and 
dotted lines show virtual levels of equidistant separation 
gfiBHo = hf, which strongly enhance the multi-photon 
transition probabilities. Since a <C hf, the distance be- 
tween the real and virtual levels in Fig. 1 is exaggerated 
for clarity. The number of arrows for each monochro- 
matic transition suggest the type of multi-photon coher- 
ent excitation (one, three or five-photon). A coherent 
Rabi superposition of the m = —5/2... 5/2 states (as 
counted on the left side of the figure) can thus be ob- 
tained, strongly dependent on the microwave drive power 
and nutation time. 



III. EXPERIMENTAL PROCEDURE 

Rabi oscillations measurements have been performed 
in a Brukcr Elexsys 680 pulse EPR spectrometer work- 
ing at about / = 9.6 GHz (A-band). The sample is 
a 3x3x1 mm 3 single crystal of MgO doped by a small 
amount of Mn 2+ (about few ppm) . The orientation be- 
tween the sample and the static field is controlled by a 
goniometer (precision 1°) with the rotation axis parallel 
to the microwave magnetic field direction h mw . In the 
measurements presented here, ft m «,||[— 110] (see Fig. 2 
insert). All measurements have been made at room tem- 
perature. Calibration of the microwave field h mw have 
been made using the BDPA standard: a small amount 
of BDPA containing isotropic S=l/2 spins give a Rabi 
frequency of exactly Fr — fiBh mw /h. 

By design, the detection is sensitive to frequency / 
and therefore it can probe only transitions between con- 
secutive spin projections, mf>m + l. As the transition 
I — 1/2) o 1 1/2) is of highest probability, as ensured by 
Fermi's golden rule, and is not sensitive to crystal strain 
effects wc choose to use the levels m = ±1/2 as a probe 
of the six-level dynamics. 
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FIG. 1. Level diagram of Mn 2+ spin in a cubic crystal field 
for a Zeeman splitting yHo — Ex/2 ~ -E-i/2- Arrows indicate 
one photon (in black), three (in red) and five (in blue) pho- 
ton monochromatic transitions, between spin projections m 
as shown on the left side. The dashed and dotted lines show 
equidistant virtual levels enhancing the multi-photon transi- 
tions. 
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FIG. 2. (Color online) Pulse sequence used in the multi- 
photon Rabi oscillation measurements, which starts with a 
strong excitation pulse Pr inducing the multi-photon nuta- 
tion, a) Right after the pulse P_r, the FID signal gives the 
transverse magnetization state between -1/2 and 1/2. b) Af- 
ter a time rp T2, a tt/2 pulse rotates the longitudinal mag- 
netization so it can be probed and so is the density of states 
in levels | — 1/2) and |l/2). (see text for more details). The 
insert shows the static and electromagnetic fields orientations 
by respect to the crystalline axes. 



By definition, the magnetization is given by (Si) = 
Tr(pSi) where p is the density matrix and the Si is the 
S=5/2 spin operator i = x,y, z. Since we probe only the 
transition | — 1/2) «-> 1 1/2) we define the magnetization of 
the subset (Si) 1/2 = T^(pi/2Si) where pi/ 2 is the central 
block matrix of p corresponding to levels ±1/2. 

The drive and detection of the multi-photon Rabi os- 
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dilations arc implemented with the pulse sequence pre- 
sented in Fig. 2. A drive pulse P_r of duration tr and 
resonant to the | — 1/2) <H> | + 1/2) transition, is applied 
at t = 0. At the end of this pulse, the density of states 
has been coherently changed. 

Right after the Pr pulse, spin dephasing and line inho- 
mogencity induce a free induction decay (FID) signal 23 . 
The first point of this FID (at t=TR_) gives the transverse 
magnetization traced on the | ± 1/2) subset: (S x )i/2 = 
P-i/2,1/2 + /°i/2,-i/2- Unfortunately, the dead time r^t of 
the spectrometer (~ 80 ns) prevents the measurement of 
this first point. If the entire inhomogenous line of width 
AHo is excited by a sufficiently short P r pulse (non selec- 
tive pulse, here t^ 1 > jAH /2ir = 0.27 MHz), the FID 
is simply the Fourier transform of the absorption line and 
is an exponential decay in MgO:Mn 2+ . Therefore, we in- 
tegrate the detected FID and the result is proportional 
to its first point. Note that for a broad line (not our 
case), as often seen in solid state paramagnetic systems, 
the FID oscillates and the first point cannot be found so 
easily 24,25 . 

The pulse sequence in Fig. 2(a) detects the coher- 
ent evolution of the transverse magnetization (S x )i/2> 
as a function of pulse length tr. To probe the coher- 
ent evolution of the longitudinal magnetization (5 f 2 ) 1 / 2 = 
Pi/2,1/2 — P-1/2,-1/2 (also traced on the | ± 1/2) subset), 
we use the sequence shown in Fig. 2(b). After the Rabi 
pulse Pr, one waits a time tp — tr, smaller than system's 
relaxation time T\ but larger than the decoherence time 
T 2 . Thus, at t = Tp, only the longitudinal magnetization 
is non-zero 26 . 

A ir/2 pulse follows, rotating the | ± 1/2) population 
in a coherent mixture, which is located in the xy plane 
(or the "equatorial plane"). The detection observes the 
free induction decay (FID) towards zero of this mixture, 
giving essential information such as: the initial total mag- 
netization in the xy plane, its FID decoherence time (T 2 *) 
and potential shift away from the resonance static field. 

As indicated with arrows in Fig. 1, in a cubic symmetry 
the three and five photon transitions do use the | ± 1/2) 
states as intermediate ones. Therefore, the ir/2 pulse 
can distinguish between the various multi-photon Rabi 
oscillations, due to their different frequencies as detailed 
in the theoretical section below. 



A. Analytical calculation 

Let us consider a quantum system with six states \S Z ), 
^={-5/2,-3/2,-1/2,1/2,3/2,5/2}, irradiated by an 
electromagnetic field in resonance with the —1/2 and 1/2 
levels. We assume initially 7 = 0, but we will describe 
the effect of the nuclear spin at the end of this section. 
The Hamiltonian of the system is : 

5/2 

H = E + V(t)= Yl E Sz \S z )(S z \+V(t), (2) 

S t =-5/2 

with Esz the static energy levels, V(t) = ^h mw (S + + 
D_) cos (27r/i), S+/S- the raising/lowering operators 
and / = Ex/2 - E-i/2- 

In a cubic symmetry and in first order in perturbation 
(a -C Ho), the static energy levels are given by: 

E ±5/2 = ±5/2 7 ff + l/2pa + 0(a 2 ) 

E ±3/2 = ±3/2 7 ff - 3/2pa + 0(a 2 ) (3) 

Kti/2 = il/27-ffo + pa + 0(a 2 ) 

where 19 p = 1 - 5(A 2 B 2 + B 2 C 2 + C 2 A 2 ), with A,B 
and C representing the cosine directors of Hq with the 
crystalline axes. For our geometry, p is given by: 

p = 1 - 5 sin 2 9 + ^ sin 4 0, (4) 

with 9 the angle between Hq and the c axis [001]. Since 
Hq h mw , the negative frequency part of V(t) respon- 
sible of the Bloch-Siegert shift 28 can be neglected and we 
can use the rotating wave approximation (RWA) to make 
Eq. (2) to be time independent. Since in a conventional 
pulsed spectrometer wc have access only to the first har- 
monic detection (hf), we define the unitary transforma- 
tion U(t) = cxp(—i2irfS z t) and apply it to Eq. (2). In 
the rotating frame, the Hamiltonian (2) becomes 21,29 : 

BJJ 

Ur W a = UHU^ + ih—U^ = (5) 



IV. RABI ROTATIONS IN A TUNABLE 
MULTI-LEVEL SYSTEM. THEORETICAL 
TREATMENT. 

The aim of this section is to present a model which de- 
scribes the coherent multiphoton Rabi oscillation that 
occur in a quasi-harmonic large spin system under 
monochromatic radiation. To compute the multiphoton 
Rabi frequencies, we used the density matrix theory 2 
applied to the Mn 2+ ion (S—5/2, 7=5/2) in the rotating 
frame approximation. 



/ \pa ^-V \ 
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where V = 7/i mM ,/2. By diagonalization, the eigenener- 
gies E n /pa of the dressed states l^n) are calculated as a 
function of h = V/pa, as illustrated in Fig. 3 (for p > 0). 

For high anharmonicity or low microwave power h — > 0, 
in which case |*i) = (|l/2) + | - l/2))/\/2 and |* 2 ) = 
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FIG. 3. (color online) Eigenvalues of Hrwa as a function of 
h. The dashed lines are the limit h mw <C p.a. 



(|l/2) — | — l/2))/-\/2 are the coherent superposition of 
states | ± 1/2). The other wavefunctions remain un- 
changed |* 3 ) = |5/2), |* 4 ) = | - 5/2), |tf 6 ) = |3/2) 
and |* 6 ) = |-3/2). 

The splitting between and l^) is detected as an 
oscillation of frequency Fr = Ei—E\ = 3^h mw , in agree- 
ment with the general formula 22 : 

F R 1] = ^h mw ^S(S + l)-S z (S z + l) (7) 

for one photon Rabi oscillations between consecutive 
states S z and S z + 1. 

When the microwave power h mw increases, the degen- 
eracy of \*f? n ) (n =3,4,5,6) is lifted, allowing coherent 
three and five photon transitions (| — 3/2) «-> |3/2) and 
I — 5/2) o 1 5/2) respectively). Note that two and four 
photon processes are out of resonance and therefore not 
discussed here. In first order perturbation in h, only the 
one photon process exists between states and | \3> 2 } - 
In third order perturbation, a three photon process mixes 
the states | ±3/2), leading to a pure 3-photon Rabi oscil- 
lation with frequency F R = ^h 3 pa, as detailed in the 
Appendix. This three photon process perturbs the afore- 
mentioned one photon process which now has the Rabi 

frequency F R ' = F R - ^h 3 pa. 

In fifth order perturbation, the five photon process 
mixes the states | ± 5/2), leading to a Rabi frequency 
F R = ^h 5 pa (pure 5-photon process). Similarly, the 
five photon process will change the dynamics of the 
three and one photon Rabi oscillations (see Appendix A, 
Eqs. A3-A5). 

For h — ¥ 00, the microwave field excites all transitions 
and the six levels are coherently driven by a 5-photon 
process. The dynamics of such a system becomes that 
of an isotropic spin S=l/2 . The diagonal of the Hamil- 
tonian (6) becomes negligible and the diagonalization 
simply rotates the Hamiltonian along the h mw field. 



When the nuclear spin is included, the model becomes 
more complicated but analytical solutions can still be 
found. The full spin Hamiltonian of 55 Mn 2+ (S = 5/2, 
7 = 5/2) is a 36 x 36 matrix. If we assume no forbidden 
nuclear transitions (Am/ = 0, with to/ quantifying the I 
projections), the Hamiltonian (1) can be separated into 
six 6x6 matrices of form Eq. (2). 

The static energies E§, now depend also on mi and 
A and have been computed in first order perturbation 
in A/ Ho by Low 19 . Having the new Eg z renormalized 
by the hyperfine field, we can use the same procedure as 
above to compute the eigenvalues of the dressed states 
and the Rabi splittings. 

Finally, note that all other transitions can be computed 
by the method presented here, by changing the value of 
/ in Eq. (2). For instance, for / = (E_ 3 / 2 — E 1 / 2 )/2 one 
can study the 2-photon Rabi oscillations in the transition 
|-3/2)^|l/2)). 



B. Numerical calculation 

The analytical model presented above describes the 
multiphoton frequencies of the Rabi oscillations very well 
but fails to describe their intensities (the time evolution 
of spin populations). In order to describe the experimen- 
tal results we developed a model using Hamiltonian (1), 
which shows the actual crystal field parameters and in- 
cludes the hyperfine coupling between the 5 = 5/2 elec- 
tronic spin and the I = 5/2 nuclear spin of 55 Mn 2+ . 

The anisotropy a and hyperfine A constants are two 
orders of magnitude smaller than the electron Zeeman 
interaction. Therefore, it is appropriate to assume that 
the orientation of the static field imposes the quantiza- 
tion axis z. In this base, Hamiltonian (1) rewrites as in 
(2), where E$ z are obtained by exact diagonalization and 
reflect the hyperfine coupling. 

Once we get the form (2) for the transformation 
in the rotating frame become obvious. Experimentally, 
the spectrometer signal represents the transient magnetic 
resonant signal in the rotating frame, obtained by mixing 
the signal reflected by the cavity (S(t) cos(2n ft)) and the 
reference arm (cos(27r/i)). This leads to a frequency- 
independent signal S(t) and a double frequency signal 
(which is filtered out). The unitary transformation is 
fixed by the mw-frequency /, with U(t) exp(— i2itS z ft) 
and using Eq. (5) one gets: 

'Hrwa = ^h mw S x + E — fS z (8) 

For a fixed mi , the static field Ho is chosen to satisfy the 
condition -E1/2 — + A = / where A represents the 

detuning of the static field, away from resonance. As an 
example, the eigenvalues of Hamiltonian (8) , for Ho 1 1 [100] 
and nuclear projection mj = 1/2, are shown in Fig. 4, 
as a function of field detuning A. In these simulations, 
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FIG. 4. (color online) Dressed states energies for i/ ||[100], 
h m w = 0.8 mT and mi — 1/2. The large anticrossing is 
the 1-photon Rabi splitting. The insets are close-ups of the 
3-photon and 5-photon Rabi splittings. 



the h mw = 0.8 mT, sufficiently strong to induce Rabi 
splitting with three and five photons, as shown in insets. 

The dynamics of Hamiltonian (8) is described by the 
time evolution of the 6x6 density matrix p(t): 



.dp 
'~dt 



[H 



RWA 



(9) 



In the rotating frame (Dirac or interaction picture), the 
solution is pit) = U p (t)poU£(t) with po the matrix density 
at thermal equilibrium : 



Po 



exp{-E/kT) 
Tr{exp(-E/kT)) 



(10) 



and U p (t) is the propagator operator: U p (t) = 
cxp(-i2n'HRWAt) ■ 

When the Rabi pulse Pr is applied, the spin population 
is coherently manipulated. At the end of the sequence a 
weak 7r/2 pulse will selectively probe the difference of 
population between states -1/2 and +1/2. The signal 
out of the spectrometer is S(t) = o-_ 1 / 2 (t) — &i/2(t)- Note 
that for a spin 5 = 1/2, S(t) = (S z )(t) 



V. EXPERIMENTAL RESULTS 

A typical MW absorbtion spectrum of MgO:Mn 2+ as a 
function of i?o, shows 6 sets (one for each mj projection) 
of 5 lines {S z -> S z + 1, S z = -5/2 . . . 3/2). In the ex- 
periments described below, the magnetic field is tuned to 
study resonances in the set corresponding to toj = 1/2. 
All experimental and theoretical results are valid to any 
other nuclear spin projection. 
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FIG. 5. (color online) Rabi oscillations and its Fourier trans- 
forms for three Ho orientations and for a microwave field 
ftmtu=1.48 mT. For Ho along the compensation angle (blue 
curves), the dynamics are the same as that of a S — 1/2 
system. For Ho 1 1 [001] (red) the microwave field is smaller 
than the effective anisotropy and only the 1-photon process 
is induced. For Ho 1 1 [1 10] (black) the microwave field and the 
effective anisotropy are comparable and the 3-photon process 
becomes visible. 



A. Tunable multiphoton Rabi oscillations by 
changing the orientation of the static field. 

The static field orientation is described by the p pa- 
rameter (see Eq. 4) and therefore can tune the reduced 
microwave field h = V/ (pa) . We measured Rabi oscilla- 
tions and corresponding Fourier transforms at constant 
microwave field (h mw =lA8 mT), for three different ori- 
entations (and p values) of Hq, as shown in Fig. 5. One 
can see a highly anisotropic behavior of the coherent dy- 
namics depending on the relative orientation between Ho 
and crystal axes. 

For i?o||[001] (p = 1, maximum value), only one fre- 
quency is observed but its value (57 MHz) is smaller than 
the one expected for the 1-photon Rabi frequency in the 
transition -1/2 o 1/2 in S=5/2 (62.5 MHz, see Eq. (7) 
and red dashed line in Fig. 5). Although the multi- 
photon Rabi oscillation is not clearly resolved, it indi- 
rectly slowing down the expected 1-photon dynamics, as 



described in Sec. IV A by the difference F 

(3) 



(1) 



= 



For ff 1 1 [110] (p = -1/4), in addition to the high fre- 
quency 1-photon process, a low frequency 3-photon Rabi 
oscillation becomes observable, due to a lower p value and 
thus a larger h. Finally, for 9 = Oca (the compensation 
angle of the cubic anisotropy: p(0ca) = 0) and hence 
h — y oo, only one Rabi frequency is observed (Fig. 5 
top). Its value is the same as the one expected for an 
isotropic spin 1/2 (dashed-dot line in fig 5). As explained 
in Sec. IV A and Appendix, a coherently driven 5-photon 
process (and only this process) is induced, leading to a 
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FIG. 6. (color online) Rabi frequency distribution of 
MgO:Mn 2+ while the static field Ho rotates in the plan (- 
110) and for a microwave field h„ lw = lA8mT. Ho=348.5 mT 
and /=9.676GHz are set to probe the transition \ms = 
-1/2, mi = 1/2 >o \m s = 1/2, mi = 1/2 >. The points 
come from the Fourier transform of experimental data. The 
error bars are the width of the Fourier transform (consequence 
of the damping). The contour plot is calculated with the nu- 
merical model given in Sec. IV B. 



S=l/2 spin dynamics. 

A detailed 8 dependence of detected Rabi frequencies, 
found by fast Fourier transformation (FFT) of coher- 
ent oscillations, is given in Fig. 6, for the same value 
of h mw . The error bars are FFT linewidths and are due 
to Rabi decay processes. The contourplot is computed 
using the numerical model described in Sec. IV B, with 
no adjustable parameters: the crystal field and hypcrfinc 
constant are from Rcf. 19 and the MW field has been 
independently calibrated using a standard. The agree- 
ment between our model and the experimental data is 
quite good. However there are some discrepancies when 
\pa\ is small (but not zero) compared to jh mw /2. In this 
case, the 4 EPR satellite lines are close to the central 
line and can overlap due to small inhomogencity caused 
by crystal field strain. Since the multi-photon Rabi fre- 
quency is highly dependent on \pa\ (specially when \pa\ 
is small), the distribution in pa induces a high distri- 
bution of multi-photon Rabi frequencies, which in turn 
overdamps the Rabi oscillation signal. 



B. Tunable multiphoton Rabi oscillations by 
changing the MW field. 

Aside the tuning method based on static field orien- 
tation, described in the previous section, an alternative 
method is based on MW drive intensity. To demonstrate 
the method, we apply a static field not far from the 
compensation angle 6qa described above, such that the 
spin system is slightly anisotropic. We chose an angle 
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FIG. 7. (color online) Rabi frequencies of Mn 2+ spin, as a 
function of drive field, obtained by FFT of detected coherent 
oscillations. The error bars are the FFT linewidths. The 
dashed line indicates the expected 1-photon dependence given 
by Eq. 7, whereas the dotted line models the isotropic S — § 
case (see Eq. A7). 



= 31 ± 1° and record coherent Rabi oscillations as a 
function of MW field h mw . The FFT peaks (Rabi fre- 
quencies) are shown in Fig. 7 by points and the error 
bars represent the FFT linewidths. The dashed blue line 
is the MW field dependence of Rabi frequency, for a pure 
1-photon transition between levels | — 1/2) <H> |l/2) of the 

5 = 5/2 system (Fr = 3%h mw , see eq. 7). The dotted 
black line is the MW field dependence of Rabi frequency 
in the case of a pure isotropic spin (Fr = ^h mw , see 
Eq. A7). 

Indeed, at low MW power, the system has the dynam- 
ics of an anisotropic S=5/2 system, as shown by the very 
good agreement between experimental points (up trian- 
gles) and the expected dependence for a 1-photon process 
(dashed line). At intermediate powers, the FFT clearly 
shows the existence of a three photon process (left tri- 
angles), as well as of a 5-photon process. The later is 
indicated by empty squares for shallow peaks in the FFT 
and by full squares, corresponding to well defined FFT 
peaks. 

A quantitative numerical description is quite difficult 
because of the extreme non linearity of h near Oca- At 
the compensation angle (6ca = 29.67°) h — > oo, but a 
small change in induces significant changes in h (e.g. 
at = 30°, h = 16.3h mw , at = 31°, h = 4.1h mw and at 

6 = 32°, h = 2Ah mw ). As the accuracy of our goniome- 
ter is 1°, a numerical quantitative description around 
9qa is not relevant. Qualitatively, one notes that for 
h m w > 0.5 mT, the one-photon Rabi frequency deviates 
from of Eq. 7, due to the formation of multi-photon 
processes. For high microwave powers (h mw > 1 mT), 
the one-photon FFT peak is barely visible (shown by 
empty up triangles) while the five photon process dom- 
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MW magnetic field (mT) 

FIG. 8. (color online) Rabi frequencies as FFT peaks (filled 
squares) of coherent oscillations for Ho along the compensa- 
tion direction of the cubic anisotropy, 8 = Oca- The contour 
plot is numerically computed, as described in Sec. IV B, while 
the dashed line shows the expected behavior for an isotropic 
S = 5/2 system (see Eq. A7). 



inates entirely the dynamics. This is indicated by the 
good agreement between the 5-photon FFT peaks (filled 
squares) and the dotted line resulting from an isotropic 
5 = 5/2 model (or a spin 1/2), as explained by Eq. A7. 

Our experiment demonstrates that a sufficiently in- 
tense microwave field can tune the spin dynamics from 
that of an anisotropic (uncquidistant) system, to that of 
an equidistant, quasi-harmonic 5-level system (equivalent 
to a two- level system) . 

C. Five-photon Rabi oscillations at the 
compensation angle. 

When the static field is applied along a direction for 
which 6 = 9q a (that is, 6 = 29.67°), the cubic anisotropy 
is geometrically compensated and therefore p = and 
h — > oo. As explained in Sec. IV A and Appendix, one 
expects a perfectly isotropic 5 = 5/2 spin system. 

We performed room temperature Rabi oscillations for 
this field orientation. The FFT of the coherent oscil- 
lations are single-peaked, as shown in Fig. 8 with filled 
squares (the FFT width is smaller than symbol's size). 
The contour plot is calculated using the numerical pro- 
cedure described in Sec. IV B, whereas the dashed line 
indicates the analytical result of Eq. A7. One observes 
an excellent agreement between the analytical and nu- 
merical models and the experimental data points (no fit 
parameter has been used). 

During a Rabi rotation for this isotropic case, the spin 
density gradually moves upward on the equidistant 6- 
level ladder and, when one reaches the highest level (S z = 
5/2), the dynamics is reversed. In contrast to a one- 
photon dynamics, here we have a collective five photon 
process, ensuring a continuous nutation within the S z = 



—5/2 ... 5/2 space. 



VI. DISCUSSION AND CONCLUSION. 

The exerimental data and the model presented here 
are using a monochromatic source to both induce and 
probe the multi-photon Rabi oscillations. Therefore, 
we restrict the types of multi-photon transitions that 
are measurable, to those which exist at the same res- 
onance static field as the one photon transition (used 
to induce and probe). In order to probe directly all 
types of multi-quanta transitions, two tunable frequen- 
cies will be needed: one to induce the multi-photon oscil- 
lations (i.e. to induce a 2-photon Rabi oscillation: f\ = 
(E 1/2 ~E_ a/2 )/2 and one to probe it f 2 = E 1/2 -E_ 1/2 ). 
Note that the experiments described here, with a single or 
two tunable frequencies, have the appeal of being appli- 
cable to on-chip studies, using magnetic detectors placed 
in strong magnetic fields 30-32 . 

In conclusion, we report on single and multi-photon 
coherent rotations in a S = 5/2 spin system featuring a 
cubic anisotropy. The anisotropy is much smaller than 
the Zeeman splittings, such as the six level scheme shows 
only a small deviation from an equidistant diagram. This 
allows us to tune the spin dynamics by cither compen- 
sating the cubic anisotropy with a precise static field ori- 
entation, or by microwave field intensity. In both cases, 
we can transit the system between an anisotropic to an 
isotropic situation, the later case showing a spin dynam- 
ics corresponding to a two-level system. The experimen- 
tal data is well explained by both a theoretical and nu- 
merical model. 
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Appendix A: Eigenvalues of dressed states 

Diagonalization of Hamiltonian (6) can be achieved an- 
alytically. Since the most relevant parameter is the value 
of the microwave field compared to the anharmonicity pa- 
rameter pa, we define the reduced field h = jh mw /(2pa). 
The eigenvalues of (6) are the solutions of the polynomial 
equations : 

V± = 6±9h+ 18ft 2 ± 15ft 3 . . . 

+ (-14 =f 12ft - 26ft 2 )£ T 12h£ 2 + 8£ 3 = GjAl) 
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FIG. 9. (color online) Rabi frequency dependence as a func- 
tion of the reduced MW-field h. Data points are numerically 
computed by solving Eq. Al, while lines are the analytical 
perturbation developments (A3)-(A5). Inset shows a close-up 
at low fields/low Rabi frequencies. 



(i) For h — > 0, in first-order perturbation, equation (Al) 
becomes : 

V± « 6 ± 9h + (-14 =f l2h)S T- ^h£ 2 + 88 3 = (A2) 

with solutions £\, 2 = 1± 1.5ft,, £3,4 = 0.5 and £5.6 = —1.5. 
The RWA wavefunctions arc now : |*i) = (|l/2) + |- 
l/2»A/2, |* 2 ) = (|l/2) - I - l/2))/V2, |* 3 > = |5/2), 



|* 4 > = I - 5/2), |* B ) = |3/2) and |* 6 ) = | - 3/2). Rabi 
oscillations occur only between states and ^2) with 
frequency: F 1 , = (£1 — £2)^ (see Eq. 7). We found 
the classical result that during the resonance, |l/2) and 
I — 1/2) are mixed by the EM- wave and the other states 
stay unperturbed. 

(ii) When h is increased, we computed the solutions of 
(Al) up to the fifth-order perturbation in h. The Rabi 
frequencies Fr = (£ + — £-)\pa\ are: 

, / ~ 24/i 3 768h 5 \ /A , 

_a f24h 3 339h 5 \ 

F « =pa {-xr-^r)> (A4) 

Fh = P a(^p). (A5) 

These equations are used to simulate the continuous 
lines in Fig. 9, which arc compared to the numerical solu- 
tions of Eq. Al (discrete data points). The inset shows a 
close-up view at low microwave powers and Rabi frequen- 
cies. One can conclude that for powers up to h ~ 1/2, 
the perturbativc analysis works very well. 

(iii) When h — > 00, Eq. (Al) becomes: 

V± ~ ±l5h 3 - 26h 2 £ =f I2h£ 2 + 8£ 3 = (A6) 
with solutions £± = h[±^; =p|; ±|]. They represent 
dressed states energies of an isotropic spin 5/2 with Rabi 
frequency: 

F 5/2 — F i/2 = -zh mw . (A7) 
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